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We study a half filled ladder of spinless fermions. We show that contrarily to a single chain, 
the ladder becomes a Mott insulator for arbitrarily small repulsive interactions. We obtain the 
full phase diagram and physical quantities such as the charge gap. We show that there is only 
a single insulating phase for repulsive interactions, regardless of the strength of the interchain 
hopping and single chain Mott gap. There is thus no confinement-deconfinement transition in 
this system but a simple crossover. We show that upon doping the system becomes a Luttinger 
liquid with a universal parameter K = 1/2 different from the one of the single chain (K = 1/4). 



I. INTRODUCTION 



One dimensional systems are one-of Jjhe few known ex- 
ample of non-fermi liquid behavioioB'El. It is thus of ut- 
most theoretical importance to understand how one can 
go from a one dimensional situation to a more conven- 
tional high (typically three) dimensional one by coupling 
one dimensional systems, allowing particles to jump from 
chain to chain. This question is far fijomJbcing elucidated 
despite several theoretical attempt&mBBtlll For com- 
mensurate one dimensional system another phenomenon 
appears: such systems are Mott insulators. This leads to 
a direct competition between interactions and hopping. 
The insulating behavior of the one dimensional system 
tends to kill the interchain hopping and thus to confine 
the electrons on individual chains. Conversely, a large 
interchain hopping destroys the one-dimensional charac- 
ter and thus weakens the Mott transition considerably, 
turning the system into a metal. This competition be- 
tween the Mott transition and interchain hopping has in 
addition to its theoretical ijnportance, implications for 



organic compoundsE£!lliijii3li3 Aat are three dimensional 
stacks of quarter filled chainscil. 

Unfortunately studying an infinite number of coupled 
chains is extremely difficult, so to understand such phe- 
nomenon it is interesting to investigate simpler systems 
with a finite number .of coupled chains. Such systems are 
the so-called ladderai3. They present the advantage to 
allow a careful study of the effec^p^ipepiiiMpbM-ieing 
tractable by powerfuL^^v^lBO-y'BEaOoy and 
numerical techniquescjO'EHlEZru. For commensurate lad- 
ders with spin, the relevance of interchain lj n ppr n £ was 
studied by renormalization group techniquesESElEij. De- 
pending on the ratio between the single chain Mott gap 
and the interchain hopping (suitably renormalized by the 
interactions) a very different flow of the single particle 
hopping was observed, reminiscent of the confinement- 
deconfinement transition expected for the infinite number 



of chains, even if in the_ladder there is no real transition 
but a simple crossoverEl In addition, interchain hopping 
was shown to drastimlLv-modify the critical properties of 
the Mott transitionE^JoEa compared to the one of a single 
chain. Despite these studies on commensurate ladders a 
detailed description of the phase diagram and of the na- 
ture of the Mott transition is still lacking. 

In the present paper we investigate these issues on a 
ladder of spinless fermions. Spinless fermions exhibit ex- 
tremely interesting behavior since a single chain needs 
a finite repulsive interaction before turning into a Mott 
insulator, contrarily to the spinful chain for which any 
repulsive interaction freezes the charge leaving only the 
spin degrees of freedom. One could thus naively think to 
be able to go from an insulating phase, dominated by the 
single chain gap, to a metallic phase even for repulsive 
interactions. In fact, quite interestingly, for the spinless 
ladder the Mott transition is pushed in the vicinity of 
the non-interacting point, invalidating this naive picture. 
Quite fortunately the fact that the Mott transition is now 
in the vicinity of the non interacting point allows to study 
it using standard renormalization group technique, and 
extract the complete properties of the transition. 

The plan of the paper is as follows. In section || we 
introdu ce t he model for the two leg spinless ladder. In 
section III wc study this model using the renormaliza- 



tion group technique. We show that the Mott transition 
occurs now for arbitrary repulsive interactions and com- 
pute the various physical parameters (charge gap, Lut- 
tinger liquid parameters) both analytically and by a nu- 
merical integration of the RG equations. We analyse the 
phase diagram in section [V. We show that in the lad- 



der the confinement-deconfinement is in fact a crossover. 
We also investigate the properties of the slightly doped 
ladder and point out the differences that exist compared 
to a doped single chain. Conclusions can be found in 
section [V|. Finally some technical details can be found in 
the appendix. 
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II. MODEL 

We start from spinless electrons on a two leg ladder, 
described by the Hamiltonian 

H = -t/SA„ a Ci+i,a + h.c.) + V }^n.i, a ni + i t a 

-tiEHi^ + h-c) (1) 

i 

where a = 1,2 is the chain index, t and t± are respec- 
tively the intra and interchain hopping, and V the repul- 
sion between nearest neighbors particles. 

To analyze the long distance properties of this model it 
is convenient. to use the boson representation of fermions 
operatorsHoa, valid in one dimension. Two basis are pos- 
sible: (i) one can start in the original chain basis and 
bosonize each chain; (ii) one can use the bonding and 
antibonding basis. Each basis has advantages and draw- 
backs and we will need both to tackle the Mott transition 
in the ladder, so we give both boson representations be- 
low. 

A Chain basis 

We refer the reader to the literature for the boson map- 
ping and recall here only the main steps to fix the nota- 
tions. Taking a linearized energy dispersion at the Fermi 
level, we use thcwfollowing expressions for right and left 
moving fermionsEJ. 



R,L 



VR,L 



(2) 



4> and II are canonically conjugate operators and II = 
\d x Q. Vcj) and V0 give respectively the long wavelength 
fluctuations of the density and current. We note wr,l 
Klein factors which one must introduce to reproduce the 
anticommutation properties of several fermion speciesu. 
a is a short distance cutoff of the order of the lattice spac- 
ing. With these operators the single chain Hamiltonian 
takes the form: 



H = H + H, 



"t" J^int "1" H u 



(3) 



H = ^-Jdx [(rfl) 2 + (d x 4>) 2 } 
= g J dx{d x (j)) 2 



H u — —2g v 



dx 



(27ra) : 



■ cos(40) 



where vf and g are respectively the bare Fermi velocity 
and the interaction. They are given by 



Vf — 2tasin(fc^a) 

(1 — cos(2kFa))aV 
9 = -5 



(4) 
(5) 



The imiklapp part, H Ul only appears in this form at half- 
fillingl3(i.c. when ikpa = 2tt) and is responsible for the 
Mott transition of a single chain. The interaction g can 
be absorbed in the quadratic part to give the Luttingcr 
Hamiltonian 



H = — I dx 

2tt 



uK(ttU) 2 + -(d^) 2 



(6) 



The parameters of the Hamiltonian are the renormalized 
fermi velocity u, the Luttinger K parameter and the non- 
universal umklapp coupling constant. For small V these 
parameters can be perturbatively computed: 



uK = vf 
u 2 

— = vf + aV—(l — cos(2fcpa)) 

K 7T 

g u = aV 



(7) 



However the description (^) is much more general and is 
valid even at large coupling provided the proper renor- 
malized coupling constant are used. In fact in the fol- 
lowing we will not assume such relations for the coupling 
constants and take free parameter (the param- 

eters may be tuned at will with for example a second 
nearest- neighbor interaction) . 

Using the bosonized expression (|J) for the single chain, 
we can write the two uncoupled chains in (uh as 



H int = g I dx{{d x 4> s f + (d x <b a f) 

dx 



(8) 



H u = ~4g u 



where we define 



(27ra) : 



■ cos 



(V8<f> a ) COs(V8(t> s ) 



2) 



The interchain hopping in (|l|) reads in this basis 
— t± — COS \/20a cos \/29 a 



(9) 
(10) 

(11) 



This basis has the advantage to treat very simply the 
umklapp term, but has the drawback not to reproduce 
easily the band picture of free fermions. 



B Two band basis 

Another basis is the bonding-antibonding band basis. 
We first diagonalize the kinetic energy in dfl) with 



1 



Ci,7T \[2 ^;2) 



(12) 
(13) 
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and the corresponding boson fields. In this basis the 
interchain hopping is diagonal: 



(14) 



however the interaction term is less simple to formulate. 
Rather than to use the bonding and antibonding boson 
fields it is again convenient to introduce the symmetric 
and antisymmetric combination that we denote now cf> p = 
{<j) a + <j)^)/\/2 and <\> a = ((f) — <j)w)/V2to distinguish them 
from the ones in the chain basis. This leads to the simple 
bosonized expression: 



t±V2 



dxd x 4> 



(15) 



The change of basis for the interaction term can be most 
easily performed using the transformation formulas for 
the total charge current . 



l/-Kd x <j> a = -1/-Kd x (j>p 
v/2 

-l/ird x (f)a = COS V2(f>cr cos V29 a 

na 



(16) 



Since the interaction terms are quadratic in currents 
we also need an operator product expansion to extract 
the most relevant operators. 



III. MOTT TRANSITION 

Let us now analyse the Mott transition. The single 
chain needs a finite strength interaction in order to be- 
come an insulator. Indeed, as can be seen from (^), for a 
single chain the umklapp term has a dimension 2 — 4K. 
It thus opens a gap in the charge sector and leads to a 
Mott insulating phase for K < 1/2 and a metallic (Lut- 
tinger liquid) phase for K > 1/2. For the t-V model this 
corresponds to V = 2t. In the ladder the presence of in- 
terchain hopping dramatically modifies this. As we will 
show, the ladder is an insulator for infinitesimal repulsive 
interactions. This can be easily seen by looking in the 
chain basis at the operators generated by t±. Obviously 
t± will generate terms such as cos(\/8^a) and cos(\^88 a ). 
Note that these terms only contain the antisymmetric 
field. However when combined with the umklapp (^) the 
cos(v / 80o) will generate a cos(v / 80 s )- This term has the 
dimension 2 — 2 K s , and in contrast with the single chain 
umklapp, is relevant for K s < 1 leading to a gap in the 
symmetric sector and hence to a Mott insulating phase. 
The perpendicular hopping thus reinforces the insulating 
character of the system. 

To go beyond this simple argument let us now inves- 
tigate the full RG flow using a two scale analysis of the 
relevant operators in the Hamiltonian. 



cosnS 2 = ~(1 — 
2 V 



-(d x (j>) 2 + cos2n</>) (17) 



The expression for the anti-symmetric current is how- 
ever misleading since it does not include Klein factors 
rcquiredJjo identify bosonic exponents to anticommuting 
fermionso. To obtain the exponential terms it is neces- 
sary to read the bosoniaed expression on the transformed 
four fermion operators^!. With these bosoii.operators the 
interaction term takes the following forrnEZI: 



H lnt = / dx(d x <f> p ) 2 



(18) 



~ J dx^nl + OM 2 ) 



Va f dx 



'-(cos 



7T- 
COS 



and the umklapp term is 

H u = „ / dxcos(\/8(/)p)(cos(v / 80 C r) +cos(V89 a )) 

(19) 

We now define the coupling constant g a associated with 
the operator cos(v / 8# (T ). 



SH = 2 9o 



dx 



(2na) 



(20) 



A Mott insulator for 1/2 < K < 1 

Let us first focus for values of K for which the single 
chain would be metallic. In that case single chain umk- 
lapp operator (||) is irrelevant. The gap in this regime 
thus results from the competition of the initial decrease 
of the umklapp constant in the flow and its subsequent 
growth after t± has blocked the transverse fluctuations. 

To investigate the gap let us write the full RG equa- 
tions, which include the umklapp term. Near the non- 
interacting point, if t± is small the flow is given by 



dfu 

dl 
dg s 

dl 
dg a 

dl 
d 9f 

dl 

dKa 

dl 

dK s 

dl 
dt± 

~dT 



(2 - 2K a - 2K s )g u 
1 



(2 - 2K s )g s 



-g u g a 



= (2 - 2K a )g a + iriK- 1 - K){t ± a) 2 + -g u g s 

= (2 - 2/K a )g f + iriK- 1 - K)(t ± a) 2 
1 



^(9f~ 9a~2g u ) 

-1 
2^r 

(2-^(K a + l/K a ))t ± 



2 -(.9 s 2 + 2<?u) 



(21) 



where we have set g/u — > g (u = 1) to keep the nota- 
tions simple. For g u = these equations reduce to the 
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ones obtained in Ref. [l7|. We have introduced couplings 
which are generated during the flow, though they are not 
present in the bare hamiltonian: 



SH 



2g a 
+2.9. 



d.v 



(27TO) 2 

dx 



; V8(f> a - 2p/ 



dx 



(27TO) 5 



•>V80 a 



(22) 



These new couplings are g a which is a density-density 
interaction between the two chains and g s which is the 
umklapp part of this interaction. The coupling gj trans- 
fers a pair of fermion from one chain to the other and 
corresponds to a Josephson coupling between the two 
chains. 

For a finite value of the renormalised ij_ that depends 
on the initial fermi velocity and interactions, one can ne- 
glect the interaction terms that couple the two bands in a 
non-resonant way. Hence we can use the above flow equa- 
tions up to a lengthscale l\ where t± reaches this finite 
value, t± ~ 0(1). Above the scale l\, it is more conve- 
nient to switch to the two band basis in order to study 
the flow. By definition of we discard any term that 
contains cos-s/8^0-, since such terms transfer momentum 
among the bands. Above l\ the flow becomes 



dK, 



1 



d! 2tt 2 ^ + 2' 9 ^ 



^ = {2-2K p -2K^)g u 



(23) 



In the o — sector the remaining cosv8^ CT term opens a 
gap. We note l a the scale when the coupling term 
cos V89 a has flowed to a value of order one, and the 
gap amplitude is evaluated as A CT = te~ l ° . Above the 
scale l a the umklapp operator becomes relevant since it 
is reduced to a simple cos - 



dg u 
dl 



(2 - 2K p )g u 



(24) 



These equations describe completely the Mott transition 
in the ladder system. An analytical solution can be given 
both in the limit of very small interchain hopping and for 
large interachain hopping but in the limit of very small 
interactions. A numerical integration of the equations 
allows to obtain the gap for arbitrary initial parameters. 



1 t± -> 

Since the umklapp is irrelevant for a single chain, one 
can replace the flow, when t± is the smallest scale in the 
problem, by 



dt ± 
~df 



dg u 
dl 
1 



(2 - AK*)g u 



(2--{K* + l/K*))t ± 



(25) 



where K* is the renormalized Luttinger parameter for a 
single <fhai». For the t — V model K* can be obtained 
exactlyc3H. When 1 1 ~ 0(1) we switch to the second 
set of flow equations (|2^). We do not know the precise 
renormalization of the umklapp constant between l\ and 
Iff . But since l a — l\ does not depend on t±, this gives 
a simple multiplicative constant. Above l a we use (|24|). 
Integration of the RG equations yields: 



,(2-4K*)h 



g u {h) = gu(0)e ( 
9u(l„) = 9u{h)e^ l "- h) 
.(Z P ) - Su(I ff )e (2 - 2r)(, '- y 



(26) 



We have used /3 as a constant that takes into account the 
variation of K a between l\ and l a as it flows to zero. This 
does not change the t± dependence. Collecting the inter- 
mediate results, we can give an asymptotic dependence 
of the total charge gap on t±: 



A p oc tj_ 



(27) 



While deriving this result we have neglected terms gen- 
erated to second order by t±, such as g a ,gf and g s . In 
the appendix ^ we show that including those terms do 
not affect the dependence of A p on t± . 



2 V -> 0, large t± 

Another interesting limit is when t± is comparable to t. 
Note that we always remain in the limit where t j_ < t in 
order to keep four points at the Fermi level, otherwise the 
problem would be the trivial one of a single filled band 
of fermions. In that case the initial flow does not exist 
and we start directly with (|23]). The umklapp operator 
has an initial dimension of 2 — AK and thus initially de- 
creases in the flow. However combined with the operator 
g u cos \f%9 a it generates from second order in the interac- 



tion expansion a term of the form g p cos v 80 p . This new 



operator is relevant, its scaling dimension being 2 - 
The flow equation of this operator reads: 



^ = (2 
dl y 



1 

2K p )g p + -^gug. 



2K, 



(28) 



Using the flow equations at large t±_ (23), we have 
g u {l) « Ve-( 2 ~ V V and g a (l) = V + 0{VW) , these ex- 
pressions are valid at the beginning of the flow and we 
have taken into account only the leading dependence on 
V. This leads to an approximate expression for g p , using 



the fact that for small interactions 2 — 2K p oc V : 



V 2 s 

>(0 = 



,-2i 



(29) 



The coupling g p starts at zero and is driven by g u which 
is decreasing rapidly. We may now determine the scale 
I* where g p becomes larger than g u . Using (|2|) one gets: 



V 2 



3 vl * = y e -(2-v)r 



(30) 
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insulator 



V^O 



t,->0 



Aich> tj_\ 



metal 



1/2 



1 K* 



FIG. 1: Phase diagram as a function of interaction and 
interchain hopping (tj_). The ladder is an insulator for 
any repulsive interaction (K* < 1), while the single chain 
is insulating for K* < 1/2 . The dashes represent a cross- 
over between a region where t± is relevant in the RG 
sense and a region where it is irrelevant (A lc ^ is the gap 
on a single chain). For 1/2 < K* < 1 the charge gap, 




FIG. 2: Typical RG flow, for 1/2 < K < 1, showing 
the two regime: below l\ (scale where t±_ becomes unity) 
and above where some couplings are cut by momentum 
conservation enabling the umklapp term to flow to strong 



A p , vanishes as a power law of t±_ for small t± and as coupling. The quantities l a and l p show the scales where 

ga/nu and g p /iru become unity at I = l a and I = l p 
respectively. The dashed curve represents the flow for 
t± = and the arrow denotes K*(~ 0.59). 



InAp oc hiV/V for small interaction (see text). 



which yields I* cx ln^, We note that at this scale g a 
has hardly changed at all, which validates the expression 
for g p . Indeed g a reaches a value of order one at l a cx y 
which for asymptotically small interactions is much larger 
than I*. The charge gap is dominated by g p and we may 
safely drop the original umklapp operator which gives 
only subdominant contributions. The gap in the total 
charge sector is given by the scale where g p reaches a 
value of order one. This gives a gap 



lnA„ oc 



InV 



(31) 



We note that the gap in the charge sector decreases faster 
than the gap in a— sector, where In A a oc V^ 1 . 



3 Numerical solution of the equations 

These asymptotic behaviors and the phase diagram are 
shown in Figure [l]. To go beyond the asymptotics either 
at small t± or at the transition V — > 0, one needs to 
numerically integrate the flow (El 
of such a flow is shown on Figure 
initial decrease of the umklapp in the initial phases of the 
flow, followed by its subsequent increase when gf ~ 0(1). 
The charge gap is given in Figure [| 





25|J24|). One example 
One clearly sees the 



FIG. 3: Charge gap as a function of t± for finite interac- 
tions characterized by the Luttinger parameter K. The 
inset shows the power law behavior of the charge gap, for 
small interchain hopping. 
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(a) 



(b) 



FIG. 4: Degenerate ground states for two uncoupled 
chains, (a) and (b). The interchain hopping selects the 
out of phase configurations (b). Such a ground state 
is the same in the whole insulating phase, regardless of 
whether A lc h > t± or not (see text). 



B Gap for K < 1/2 

The picture is quite different in this case since with- 
out interchain hopping a gap would be present on each 
chain. For small t± the model is conveniently studied in 
the chain basis, using the flow j2l|). Because the field <pi t 2 
orders, the single particle hopping between the chains is 
irrelevant (for small t±). The only remaining relevant 
coupling is thus the generated interchain density-density 
interaction g a . This coupling opens a gap in the antisym- 
metric sector and orders <p a . The physics is quite clear: 
in the absence of interchain coupling each chain has a 
charge density wave ground state that is double degener- 
ate. In the presence of t± , particles still cannot hop from 
chain to chain because of the Mott gap in each chain, 
but the virtual hops tend to lock the two CDW relative 
to each other. This is shown on Figure ||. The Mott gap 
is thus essentially here the single chain gap. The gap as 
a function of K obtained by numerical integration of the 
flow is shown in Figure |^. 



IV. PHASE DIAGRAM 



The results of section [II allow to draw the phase dia- 
gram of the commensurate ladder, as shown on Figure [l]. 
The very existence of an insulating phase in the ladder 
for 1/2 < K < 1 prompts for several questions. 



A Confinement vs crossover 

Apparently we have to face two very different behaviors 
depending on the strength of t±. For K < 1/2 and smal l 
t± , t± renormalizes to zero as was shown in section III B , 



whereas for large t± (or for 1/2 < K < 1 where the single 
chain gap is absent) t±_ renormalizes to large values. One 
thus seems to have a confinement-deconfinement transi- 
tion, similar to the one expected for an infinite number 
of chains, induced by the competition between t± and 
the single chajnJVIott gap Ai c h. The change of behavior 
occurs whencJ'Ej 
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FIG. 5: Charge gap for finite t± as a function of the 
Luttinger parameter K, the inset shows curves for differ- 
ent values of t±. For K < 1/2 the Mott gap is essentially 
the single chain gap, whereas it is strongly t± dependent 
for 1/2 < K < 1. 



where t± is the renormalized interchain hopping. We 
thus have to determine here whether we have two dif- 
ferent insulating phases. This can be done by looking 
at the fields that order in e ach ca se. In the "confined" 
phase we showed in section III B that both (j) s and 4> a 
acquire mean val ues. I n the "deconfined" phase, as was 

and <j)p orders. The physical 



shown in section III A 



observable corresponding to an out of phase charge den- 
sity wave takes the following form with the two different 
set of boson operators: 



J. 



(33) 



= — sin \/2(j) s 



sin 



V2<j> a 



It 



is easy 

COSV^ct) = 
COsV8</> a ) = 



= — sin \/2(f>„ sin V29 a 

TTCL 

to see from (|3|) that the conditions 
— 1, (cos V8(j)p) = — 1 on the one hand and 
: — 1, (cos V^s) = — 1 on the other both 
give a non zero value to (J| & ). Thus the two "phases" 
have the same ordered fields (written in a different ba- 
sis) and there is no transition but a simple crossover. 
The out of phase charge density waves ground state that 
arises naturally from the K < 1/2 picture and is shown 
in Figure || stays valid even in the other limit. 



t 



cir 



Vlch 



(32) 
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TABLE I: Luttinger liquid parameters at the Mott 
transition, when varying the strength of the interaction 
(Mott-U) and when varying the doping (Mott-5). 



Parameter 


lchain 


Ladder 


-^Mott-U 


1/2 


1 


"Mott-<5 


1/4 


1/2 



B Doped ladder 

Let us now determine the phase diagram of the ladder 
system away from half filling. Introducing a chemical 
potential would change the Hamiltonian into 

H = H 1/2ming - V2fi/n J dxVc^p (34) 

As we have seen that the operator opening the charge 
gap is always 

#Mottladder = J dx COs(V8,(j) p ) (35) 

instead of 

HMott lchain = J dxcos(4:(f)p) (36) 

for the single chain, tteing the known generic solution 
of such Mott systemsEI] we can see that in the ladder 
the doping will destroy the Mott phase when fi = A. 
The properties at the transition are the generic proper- 
ties of the Mott transition in one dimension when varying 
the doping (Mott-<5), namely: (i) a dynamical exponent 
z = 2; (ii) a divergent compressibility on the metallic 
side. The luttinger liquid exponent is universal and is 
K* = 1/2. This values is quite different from the one 
of a single chain. The differences are recalled in Table |. 
The difference in K* should be observable in numerical 
calculations such as exact diagonalization and dmrg. It 
might be a more clear signature of the difference between 
the ladder and the single chain that the gap itself spe- 
cially for K < 1/2. 

These different values of K* correspond to two differ- 
ent kinds of elementary excitations. For the gapped sin- 
gle spinless chain the elementary excitations are domain 
walls separating charge density waves with a different 
phase, the equivalent of the excitations called spinons 
in the spin— 1/2 Heisenberg model. In a path-integral 
picture these corresponds to configurations where the 
field <p is step-like with for instance <j>{— oo, t) = and 
</>(+oo,<) = 7r/2. However such a configuration on one 
chain with a nearly constant configuration on the sec- 
ond chain has an infinite cost when there is a cos \f8<fr p 
term in the model, no matter how small the coupling 
constant. Therefore the low-energy excitations allowed 
for two coupled chains are kinks (two domain walls con- 
fined in the same region). This behavior is reminiscent of 




FIG. 6: Phase diagram as a function of the chemical 
potential /i. The dashed line shows the metal-insulator 
(commensurate-incommensurate) transition for a single 
chain (beginning at finite interaction V — 2t), the con- 
tinuous line corresponds to the spinless fermionic ladder. 
K* is the universal value of the luttinger liquid parameter 
at the transition. 

the double-sine-Gordon model that appears for a spon- 
taneously dimcrized frustrated spin— 1/2 chain (next- 
nearest-neighbor exchange) with an infinitesimal bond 
alternation perturbationEH. 

A summary of the phase diagram is shown on Figure ra. 



V. CONCLUSION 

We have investigated in this paper commensurate lad- 
ders of spinless fermions. We have shown that contrarily 
to the single chain, that needs a finite repulsion to be- 
come a Mott insulator, a two leg ladder turns insulating 
for arbitrary small repulsion between the fermions. The 
interchain hopping thus paradoxically reinforces the in- 
sulating behavior of the system. We have computed the 
charge gap as a function of the Luttinger exponent K of 
a single chain. For values of K for which the single chain 
system would be metallic (1/2 < K < 1) the charge gap 
vanishes as a power law of the interchain hopping ij_ for 
small t±, and faster than the standard BKT behavior, 
usually characterizing metal insulator transitions in one 
dimension, when the interactions tend to zero. The inter- 
chain hopping also affects the universal properties of the 
transition that occur upon doping the system. Indeed for 
very small doping the Luttinger liquid parameter takes 
the universal value K* — 1/2 instead of K* = 1/4 for a 
single chain (when a Mott gap is present). This differ- 
ence in Luttinger liquid parameters should be accessible 
in numerical simulations such as exact-diagonalization or 
dmrg. 
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Given the fact that when a Mott gap is well estab- 
lished in a single chain the interchain hopping scales to 
zero, whereas it scales to large values otherwise, it was 
important to know whether there was a confincmcnt- 
deconfincment transition. We have shown that there is 
in the spinless ladder only a crossover and that the two 
seemingly different insulating phases (large Mott gap on 
a single chain and small t± and large t± and small (or 
zero) single chain Mott gap) are in fact of the same na- 
ture. 



This leaves open how this cross-over should evolve tor 
an increasing number ol coupled chains to give back the 
confinement-deconfinement transition expected for an in- 
finite number of chains. 
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APPENDIX A: TERMS GENERATED TO 
SECOND ORDER BY T± 



The RG equations ( |21| ) show that during the first part 
of the flow up to li, t± generates to second order the term 
g a cos y8(j) a . Contracting g a with g u gives a contribution 
to the renormalization of the term g s cos s/8(/) s (noticing 
that <f) p — 4> s , we switch to the notation g p instead of 
g s ). This term is initially not present in the hamiltonian. 
However it is generated and relevant for any repulsive 
interaction. One must thus make sure that this term does 
not grow so fast as to open the gap before the original 
umklapp term. If it were the case this would change the 
expression for the charge gap compared to the one given 
in the text. 

To determine the growth of this coupling, we have to 
determine g a (l)- Naively the scaling dimension of g a is 
2 — 2K, however it is generated by t\, term of dimension 
4 — (K + 1/K). Comparing the two rates shows that 
the source term dominates. Integration then shows that 
g a (l) oc tj^Q). We may now repeat the same analysis 
for g p , indeed its scaling dimension is 2 — 2K but it is 
actually driven by g u {l)ga(l), term of dimension 6 — 5K — 
1/K, dominant in the regime considered here (K > 1/2). 
Consequently one has g p (l) oc g u (l)g a (i)- Thus at the 
scale li, where t± is of order one, one has g P {h) ~ 9u(h)- 

To complete the argument we note that the couplings 



9 P (l) cos 



and g u (l) cos 



cos 



are above l\ 



basically the same operator. Indeed in the a— sector g a 
has flowed to strong coupling above l a , (cos-\/8# CT ) ^ 0. 



Recalling that the flow between li and l a does not change 
the asymptotic dependence on t± and that g u (h) an d 
g p (h) are proportional we recover by following the flow 
of the coupling g p the behavior explained in the text. 
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